Diffusion and transfer of entanglement in an array of inductively coupled flux qubits 
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A theoretical scheme to generate multipartite entangled states in a Josephson planar-designed 
architecture is reported. This scheme improves the one published in [Phys. Rev. B 74, 104503 
(2006)] since it speeds up the generation of W entangled states in an M x N array of inductively 
coupled Josephson flux qubits by reducing the number of necessary steps. In addition, the same 
protocol is shown to be able to transfer the W state from one row to the other. 
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I. INTRODUCTION 

In the last years, condensed- matter architectures 
based on Josephson-junction qubits have appeared 
to be promising candidates for quantum information 
processors^ These solid-state systems can be scaled up 
to a large number of qubits and satisfy DiVincenzo's re- 
quirements for quantum computing,— i.e. state prepara- 
tion, manipulation, and readout. Among them, qubits 
on the basis of the superconducting quantum interfer- 
ence devices (SQUID) offer the possibility of realizing 
switchable (inductive) interbit couplings 4 ^^ and there- 
fore provide promising platforms to generate multipartite 
entanglements among "macroscopic" quantum systems 
in deterministic ways. Remarkable experimental achieve- 
ments with flux qubits include the realization of complex 
singlc-qubit operation schemes^ the generation of entan- 
gled states^ in systems of coupled flux qubitSfi^ while 
the direct measurement of the entanglement in Josephson 
architectures has been performed via state tomography-ii 
for two superconducting phase qubits. 

Within such frameworks, we proposed a scheme for the 
generation of a W entangled state in a chain of N spa- 
tially separated flux qubits by exploiting their sequential 
couplings with one of them playing the role of an en- 
tanglement mediator^ We remind that the TV-partite 
W entangled state is a natural generalization to TV-qubit 
state of the tripartite W state \W) 3 = (jlOO) + |010) + 
|001))/>/3, i.e., the totally symmetric (apart from possi- 
ble phase factors) quantum superposition of N two-state 
systems where only one of them is in its excited state. 
In Ref. HH, the success of this scheme relies on the pos- 
sibilities of both preparing the initial state of the qubits 
and tuning the coupling energy and/or the interaction 
time between each qubit and the mediator, provided the 
time necessary for the desired quantum processes is short 
enough with respect to the decoherence time. 

The protection against noise is evidently one of the 
central issues in quantum information technology and the 



reduction of the duration spent for specific quantum op- 
erations is important for it. In this paper, we improve the 
scheme proposed in Ref. [ID by analyzing the dynamics of 
an array of flux qubits, which can be selectively coupled 
in pairs for instance by exploiting the tunable flux trans- 
former proposed by Castellano et al.fi and demonstrate 
how it is possible to diffuse a W state prepared in one 
row to two or more rows with a few steps. Such a scheme 
helps to reduce the time for the generation of multipartite 
entanglement. Furthermore, we show that the same pro- 
tocol also provides us with a way to shift, or transfer, the 
W state from row to row. We emphasize that W states 
are promising candidates for the experimental realization 
of quantum information processing in multiparticle sys- 
tems since they possess entanglement robustness against 
local operation even under qubit loss. 



II. AN ARRAY OF QUBITS 

The idea presented in this paper is based on the theo- 
retical proposal in Ref. [H to generate a W state in a chain 
of rf-SQUID or persistent current (3JJ) qubits. We first 
recapitulate its essential idea within the present setup il- 
lustrated in Fig.[TJ Here, in order to minimize the suscep- 
tibility to external noise of a large-inductance rf-SQUID, 
as proposed by Mooij et aL,— we consider the planar 
array sketched in Fig. 1 constituted by M x N spatially 
separated, and consequently not directly interacting, tun- 
able 3JJ qubits, that is a superconducting loop contain- 
ing three Josephson junctions, two of equal size (i.e. with 
Ej t i = Ej t 2 = Ej) and the third one smaller by a factor 
a (i.e. with Ej$ = aEj, a < 1). This parameter may 
be adjusted for instance by substituting the third junc- 
tion with a dc-SQUID behaving as an effective JJ with 
tunable Josephson energy, Ej^ = Ej^(4> c ), <f> c being an 
additional control flux threading the dc-SQUID loop. In 
such conditions, by applying an external flux <j) x close 
to a half-integer number of flux quanta, 4>q = h/2e, and 
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FIG. 1: Schematic illustration of an M x TV array of Joseph- 
son flux qubits. An inductive qubit-qubit coupling is realized 
by means of a superconducting switch, namely a transformer 
with variable flux-transfer function 7?.(</)c™ Tl '), as proposed in 
Ref.|. It is possible to control the flux-transfer function, and 
therefore the inductive-coupling constant, by modulating the 
critical current of the inner dc-SQUID of each transformer via 
an externally applied magnetic flux $^ n \ Each individual 
coupling between a pair of qubits is effectively turned on by 
adjusting the control fluxes of the relevant "switches" with all 
the other "switches" kept off. 



choosing a as 0.8 the potential energy of the total system 
forms a double well which permits two stable configura- 
tions of minimum energy corresponding to two persistent 
currents ±I p « ±2naEj/(j)Q in the loop. This fact allows 
to engineer a two-state quantum system (qubit) whose ef- 
fective Hamiltonian, in the basis of the two energy eigen- 
states |0 mn ) and |l mn ) of the (m,n) qubit of the array 
(which at <j) x = 4>q/2, are maximal superpositions of the 
two persistent-current states \L mn ) and \R mn )), reads 



H mn = -fuv mn a^ mn \ uj mn = yj A^„(</j c ) + e mn ((j> x ), 



(1) 

where m = 1, . . . , M, n = 1,...,N, and ai™ 1 ™^ = 
|lmn)(lmn| - |0„ m )(0„ m | is a Pauli operator for the 
(m, n)th qubit. The energy spacing, corresponding to 
a transition frequency ui mn typically in the range of mi- 
crowaves, can be tuned by properly selecting both the 



tunneling frequency, A 



mn , between \L mn ) 



and \Rmri) and 



b x ) = 2I P 



6o/2)/ft both depending on the sys- 



tern parameters. 

In the following discussion, we assume that all the 
qubits have a common energy gap u> = u) mn , V(m, n). 

The qubits are coupled with each other as depicted 
in Fig. [1] The coupled dynamics of the total system is 
described by the Hamiltonian 



H = Hn 
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where 
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mn.kl\' J + a - 
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(3) 



is the rotating-wave coupling between qubits (m, n) and 



(k,£) with a y l 



(mn) 



|l m n)(0 m „| and a 



(ran) 



| Omn) (It 



the raising and lowering operators for qubit (m, n), cor- 
respondingly. 

The coupling constants gn,in = 2(^n j i„/L)(0o/2) 
(n = 2, . . . , TV) between the first and the nth qubit in 
the first row can be turned on and off via controlling the 
magnetic fluxes 4>cx externally applied to the (m, n) 
qubit as shown in Fig.[T] Analogously we may control the 
coupling constants g mn ,( m +i)n = 2(& n „,( m+ i)„/£)(>o/ 2 ) 
(with m = 1, . . . , M — 1, n = 1, . . . , TV) between the 
mth and the (m + l)th qubits in the nth column. 

Here (u,i„ = r n K(^')Tl((t)^)r ln and £ m „,( TO+ i)„ = 

^(m+i)n? ^*mn being the flux transform- 
ing ratio between the arm of the transformer and the 
qubit (m,n). In the following calculations, we work in 
the interaction picture with respect to Hq. 

If the inter-row couplings are turned off, q m „ („, \\ 7I = 
0, the system is essentially the one analyzed in Ref.Qll It 
is therefore possible to generate an TV-partite W state in 
the first row, among qubits (1, 1), . . . , (1, TV), as follows: 

• We prepare the initial state |Wo) = |ln • • • Oijv) <8 
|0 2 i • ■ • 2 n) ® • • • <S> \0mi ■ ■ ■ Omn), with only the 
qubit (1,1) in the excited state and the rest of the 
qubits of the array in their own ground state. 

• The coupling 511.12 is turned on during a proper 
time interval < t < t\ with other couplings off. 

• The coupling 311,12 is turned off at t = T\ and in- 
stead 311,13 is on during t\ < t < t\ + t%. 



• In this way, qubit (1,1) 
(1,2), . . . , (1,TV) one by one. 



is coupled with 



It is shown in Ref. [l2j that, by exploiting the knowledge 
of the coupling constants gxiu n +i) when the relevant 
interaction is turned on, by setting the interaction times 
t„ so that 
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cos( 



TV-n 



(4) 



y/N-n+l VTV-n+1 
with 9 n = gii,i(n+i) T n/fc (n = 1, . . . , TV — 1), a W state 
K) = ® IO21O22 ■ • ■ 2N ) ® • ■ 
with 

\W) m = -^= (|l m iO m2 . . -0 mN ) - i|0 ml l m2 • • .0 
V TV \ 

— • • • — «|O m lO m 2 . . . lm«)j 



|0mi0m2 • ■ ■ Omn) 
(5) 



1 N 



(6) 



is generated in the first row m — 1 of the array. 
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We underline that rapid-singlc-flux-quantum (RSFQ) 
Josephson-junction based logic circuit al 15 make it possi- 
ble to produce flux pulses characterized by rise/fall times 
t r I f of the order of 10 ps. Therefore we are able to obtain 
switching times much smaller than the duration of any 
step in our scheme, typically less than or of the order of 
the inverse of coupling energy h/gj « 2 ns. 

III. DIFFUSION AND TRANSFER OF W STATE 

Let us now discuss schemes for spreading the W state 
([5]) prepared in the first row of the array to other rows, as 
well as for transferring it from row to row, by making use 
of the switchable inductive coupling g mn ,(m+i)n between 
each qubit of the mth row (m = 1, . . . , M — 1) and the 
corresponding one in the (m + l)th row. To this end, we 
consider the following "collective" step by step scheme: 

• Each of the TV qubits in the first row, already pre- 
pared in the W state (0), is put in inductive in- 
teraction with the corresponding one in the second 
row during < t < Ti, by turning on the couplings 
5n,2i = .912,22 = • • • = 9in,2N = ,9i ^ 0, while 
other qubits evolve freely. 

• Ati = ti, the couplings gii.2i, 312,22, ■ • - ,9in,2N are 
turned off, and 321,31 = .922,32 = • • • = g 2 N,m = 
3n 7^ are turned on for time interval t\ < t < 
t~\ + Tij, hi order to couple the qubits in the second 
row with the corresponding ones in the third. 



• Similarly, the interactions between adjacent rows 
are successively switched on and off. 

Then, by properly selecting the interaction times (by 
turning on and off the coupling constants), it is possi- 
ble to transfer the W state from row to row or to diffuse 
it to multi-rows. To illustrate the mechanism, let us look 
at a 3 x 3 array (M = N = 3) for the sake of simplicity. 
In this case, the W state prepared in the first row, Eq. 
©, reads 

hi) = \W)i ® IO21O22O23) <8> |0 3 i0 3 2033)- (7) 

By switching on the inductive couplings 311.21 = 312,22 = 
513,23 = 3i 7^ for a time period n, the W state \w\) in 
J7|) is driven into 

|*i) = (cos(9i|1ii0i 2 0i 3 ) <g) IO21O22O23) 

- icos6»i|0iili 2 0i 3 ) ® IO21O22O23) 
-icosft|0 u 0i 2 li 3 ) ® IO21O22O23) 

- isin0i|OiiOi 2 Oi3) ® II21O22O23) 

- sin0i|O u Oi 2 Oi 3 ) <g> IO21I22O23) 

- sin0i|O u Oi 2 Oi3) ® IO21O22I23)) 

® IO31O32O33), (8) 

and this is further converted into 



|*i) = (cos0i|1iiOi 2 Oi 3 ) ® IO21O22O23) ® IO31O32O33) - icosf?i|0iili 2 0i 3 ) <8> IO21O22O23) ® IO31O32O33) 

- icos6»i|0ii0i 2 li 3 ) ® IO21O22O23) ® IO31O32O33) - isineicosftilOnO^Oia) ® II21O22O23) ® IO31O32O33) 
-sin^cosflnlOnOiaOia) ® |0 2 il 22 2 3) ® |0 31 0320 3 3) - sin^ cos6>ii|0 u 0i 2 0i 3 ) ® |0 21 22 l2 3 ) ® IO31O32O33) 
-sin^ sin^nlOnOiaOia) ® |0 2 i0 22 2 3) ® II31O32O33) + i sin ft sinfti|0 a 0i 2 0i 3 ) ® IO21O22O23) ® IO31I32O33) 

+ isinft sinfti|0ii0i 2 0i 3 ) ® IO21O22O23) ® |0 3 i0 32 l3 3 )) 

(9) 

I 



after the second step with the couplings 3 2 i.3i = 3 22 .32 = 
523,33 = 3ii 7^ turned on while others are off, where 
ft/11 = 5i/ii r i/n/fr- 

Equation (|8]) clearly shows that the tuning sin ft = 

cos ft = 1/ -s/2 realizes a one-step diffusion of the W state 
up to the second row, namely from the tri-partite to a 
hexa-partite W state 

|W 2 ) = (|ln0i 2 0i 3 ) ® IO21O22O23) 
- i|0iili 2 0i 3 > ® IO21O22O23) 
-i|0n0i 2 li 3 ) ® |0 2 i0 22 23 ) 



-i|0n0i 2 0i 3 ) ® |l2i0 22 2 3) 

- |0n0 12 13 ) ® IO21I22O23) 

- I0n012013) ® 1021022123)) ® 1031032033)- 

(10) 

Or alternatively, if we select a different interaction 
strength so as to satisfy sin ft = 1 , another W state 

\w 2 ) = -i|0 u 0i 2 0i 3 ) ® \W) 2 ® IO31O32O33) (11) 

is established, that is, the tri-partite W state in the first 
row is shifted, or transferred, to the second row after the 
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one step. 

The second step further diffuses or transfers the W 
state to the third row. Indeed. Eq. ([9]) shows that 
one choice sin#i = a/2/3, cos 9\ = 1/V3 and sin#n = 
cos 0n = l/y/2 yields a complete W state all over the 
3x3 array 
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namely, the W state (© is diffused from the three qubits 
to the nine ones with only two steps; more generally, the 
W state prepared in the first row is diffused all over the 
M x N array via N — I steps. Another choice sin#i = 
sin On = 1 generates 

|tua) = -IO11O12O13) ® IO21O22O23) ® \W) 3 (13) 
transferring the W state ||SJ) in the first row to the third. 

IV. CONCLUSIONS 

In summary, we have extended the theoretical scheme 
for the generation of a W entangled state proposed in 
Ref. [l3 to a scheme for an M x N array of qubits. A re- 
markable feature is that an entanglement realized as a W 
state in one row of the array is diffused to two rows and 
the number of qubits involved in the W state is doubled 
after a single step. The entanglement is further diffused 
all over the array by repeating similar processes to yield 
an (M x iV)-qubit W state. This procedure would facili- 
tate the generation of a large-scale multipartite W state 
with fewer steps, and as a result, would help to save time 
for its generation. Furthermore, we have demonstrated 
the possibility of transferring the W state prepared in one 
row to another at will. We have illustrated these schemes 
in the context of the inductively coupled flux qubits, but 
they are also applicable to other systems, provided the 
qubit-qubit couplings are controllable. 



The experimental realization of our proposal is pos- 
sible, although it has to face various technological chal- 
lenges. The fabrication of the proposed circuit is not diffi- 
cult. Moreover, to exploit tunable qubits and flux trans- 
formers allows both to strongly reduce the single-qubit 
parameter disorder and at the same time to select com- 
mon energy gaps with a disorder of a few percent^ and to 
control the duration of intcrqubit couplings. In addition, 
state preparation can be accurately realized with well- 
defined procedures including relaxation and single qubit 
rotations. A delicate point concerns the system time 
scales against undesired effects due to the coupling of 
each qubit with bosonic baths, traceable back to the pres- 
ence of many unavoidable noise sources. For instance, the 
effective impedance characterizing the dissipative elec- 
tronic circuitry coupled to the single qubit progressively 
degrades its coherent evolution. The consequent impact 
on the single qubit decoherence rates and on the perfor- 
mance of a gate of two inductively coupled qubit have 
been studied bringing to light that, by carefully engi- 
neering the environmental impedances, the bipartite sys- 
tems is characterized by rates in the range 10 -7 -10 -6 
Since the observed relaxation and decoherence times for 
a single flux qubit are in the range 1-10 /ns2 ' 15 i 18 the pas- 
sage to a bipartite system speeds up the decoherence pro- 
cess. Thus the present scheme, shortening significantly 
the generation time of a multipartite W state, provides 
an effective way to anticipate the occurrence of coher- 
ence loss effects. We assume that the eigcnfrcqucncy u 
of a Josephson qubit is of the order of 10 GHz and that 
the inductive qubit-qubit coupling constant is of the or- 
der of 0.5 GHz<2£ Under these conditions, the length of 
a generic step (during which only a fraction of a Rabi os- 
cillation takes place) is of the order of 2 ns. This means 
that the W state (fT2")) for the 3x3 array is generated ap- 
proximately after 8 ns, which is much shorter than 16 ns 
required for the generation of a 9-partite W state by the 
scheme discussed in Ref. [I2T Extending this argument to 
a larger array (for instance, to a 5 x 5 array), we find 
that the generation of an entangled state of 25 qubits 
(which requires 8 steps) is roughly compatible with the 
currently observed relaxation times, ks 300 ns, character- 
izing two coupled qubits. We wish to emphasize that 
this estimation is suitable in our case, since during each 
step we deal with non interacting bipartite systems of 
inductively coupled qubits. If we access only one qubit 
at a step to entangle multiple qubits, such a large-scale 
entanglement cannot be established within such limited 
time. 
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